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Abstract 

Motivated by recent models proposing a hidden sector with ~ GeV scale force carriers, we present 
a search for a narrow dilepton resonance in 4 lepton final states using 536 fb -1 collected by the 
BABAR detector. We search for the reaction, e + e~ — ► WW — > where the leptons 

carry the full 4-momentum and the two dilepton pair invariant masses are equal. We do not observe 
a significant signal and we set 90% upper limits of a(e + e~ — ► WW — > e + e~e + e~) < (15 — 70) ab, 
a(e + e~ -» WW -> e+e-/z + /0 < (15 - 40) ab, and a{e + e~ -» WW -► fA + fTfJ, + fT) < (11 - 
17) ab in the W mass range between 0.24 and 5.3 GeV/ c 2 . Under the assumption that the W 
coupling to electrons and muons is the same, we obtain a combined upper limit of <r(e + e~ — > 
WW — > l + l~l ,+ l'~) < (25 — 60) ab. Using these limits, we constrain the product of the SM-dark 
sector mixing and the dark coupling constant in the case of a non-Abelian Higgsed dark sector. 



Submitted to the XXIV International Symposium on Lepton Photon 
Interactions at High Energies, 
17 August — 22 August 2009, Hamburg, Germany. 



SLAC National Accelerator Laboratory, Stanford University, Stanford, CA 94309 
Work supported in part by Department of Energy contract DE-AC02-76SF00515. 



The BABAR Collaboration, 

B. Aubert, Y. Karyotakis, J. P. Lees, V. Poireau, E. Prencipe, X. Prudent, V. Tisserand 

Laboratoire d'Annecy-le-Vieux de Physique des Particules (LAPP), Universite de Savoie, CNRS/IN2P3, 

F- 74941 Annecy-Le-Vieux, France 

J. Garra Tico, E. Grauges 
Universitat de Barcelona, Facultat de Fisica, Departament ECM, E-08028 Barcelona, Spain 

M. Martinelli afc , A. Palano ab , M. Pappagallo afc 
INFN Sezione di Bari a ; Dipartimento di Fisica, Universita di Bari b , F70126 Bari, Italy 

G. Eigen, B. Stugu, L. Sun 
University of Bergen, Institute of Physics, N-5007 Bergen, Norway 

M. Battaglia, D. N. Brown, B. Hooberman, L. T. Kerth, Yu. G. Kolomcnsky, G. Lynch, I. L. Osipcnkov, 

K. Tackmann, T. Tanabe 

Lawrence Berkeley National Laboratory and University of California, Berkeley, California 94720, USA 

C. M. Hawkes, N. Soni, A. T. Watson 
University of Birmingham, Birmingham, B15 2TT, United Kingdom 

H. Koch, T. Schroeder 
Ruhr Universitat Bochum, Institut fur Experimentalphysik 1, D-44780 Bochum, Germany 

D. J. Asgeirsson, C. Hearty, T. S. Mattison, J. A. McKenna 
University of British Columbia, Vancouver, British Columbia, Canada V6T 1Z1 

M. Barrett, A. Khan, A. Randle-Conde 
Brunei University, Uxbridge, Middlesex UB8 3PH, United Kingdom 

V. E. Blinov, A. D. Bukin0 A. R. Buzykaev, V. P. Druzhinin, V. B. Golubcv, A. P. Onuchin, 
S. I. Serednyakov, Yu. I. Skovpen, E. P. Solodov, K. Yu. Todyshev 

Budker Institute of Nuclear Physics, Novosibirsk 630090, Russia 

M. Bondioli, S. Curry, I. Eschrich, D. Kirkby, A. J. Lankford, P. Lund, M. Mandelkern, E. C. Martin, 

D. P. Stoker 

University of California at Irvine, Irvine, California 92697, USA 

H. Atmacan, J. W. Gary, F. Liu, O. Long, G. M. Vitug, Z. Yasin 
University of California at Riverside, Riverside, California 92521, USA 

V. Sharma 

University of California at San Diego, La Jolla, California 92093, USA 

C. Campagnari, T. M. Hong, D. Kovalskyi, M. A. Mazur, J. D. Richman 
University of California at Santa Barbara, Santa Barbara, California 93106, USA 

1 Deceased 



2 



T. W. Beck, A. M. Eisner, C. A. Heusch, J. Kroseberg, W. S. Lockman, A. J. Martinez, T. Schalk, 
B. A. Schumm, A. Seiden, L. Wang, L. O. Winstrom 

University of California at Santa Cruz, Institute for Particle Physics, Santa Cruz, California 95064, USA 

C. H. Cheng, D. A. Doll, B. Echenard, F. Fang, D. G. Hitlin, I. Narsky, P. Ongmongkolkul, T. Piatenko, 

F. C. Porter 

California Institute of Technology, Pasadena, California 91125, USA 

R. Andreassen, G. Mancinelli, B. T. Meadows, K. Mishra, M. D. Sokoloff 

University of Cincinnati, Cincinnati, Ohio 45221, USA 

P. C. Bloom, W. T. Ford, A. Gaz, J. F. Hirschauer, M. Nagel, U. Nauenberg, J. G. Smith, S. R. Wagner 

University of Colorado, Boulder, Colorado 80309, USA 

R. Ayad@ W. H. Toki 

Colorado State University, Fort Collins, Colorado 80523, USA 

E. Feltresi, A. Hauke, H. Jasper, T. M. Karbach, J. Merkel, A. Petzold, B. Spaan, K. Wacker 
Technische Universitdt Dortmund, Fakultat Physik, D-44221 Dortmund, Germany 

M. J. Kobel, R. Nogowski, K. R. Schubert, R. Schwierz 
Technische Universitdt Dresden, Institut fur Kern- und Teilchenphysik, D-01062 Dresden, Germany , 

D. Bernard, E. Latour, M. Verderi 
Laboratoire Leprince-Ringuet, CNRS/IN2P3, Ecole Poly technique, F-91128 Palaiseau, France 

P. J. Clark, S. Playfer, J. E. Watson 
University of Edinburgh, Edinburgh EH9 3JZ, United Kingdom 

M. Andreotti Qb , D. Bettoni a , C. Bozzi a , R. Calabrese ab , A. Cecchi afc , G. Cibinetto ab , E. Fioravanti afc 
P. Franchini ab , E. Luppi ab , M. Munerato afc M. Negrini afc , A. Petrella ab , L. Piemontese Q , V. Santoro ab 

INFN Sezione di Ferrara a ; Dipartimento di Fisica, Universita di Ferrara b , 1-44100 Ferrara, Italy 

R. Baldini-Ferroli, A. Calcaterra, R. de Sangro, G. Finocchiaro, S. Pacetti, P. Patteri, I. M. Peruzzijf] 

M. Piccolo, M. Rama, A. Zallo 

INFN Laboratori Nazionali di Frascati, 1-00044 Frascati, Italy 

R. Contri Qb , E. Guido ab , M. Lo Vetere afc , M. R. Monge ab , S. Passaggio a , C. Patrignani ab , E. Robutti Q , 

S. Tosi Qb 

INFN Sezione di Genova a ; Dipartimento di Fisica, Universita di Genova b , 1-16146 Genova, Italy 

M. Morii 

Harvard University, Cambridge, Massachusetts 02138, USA 

A. Adametz, J. Marks, S. Schenk, U. Uwer 
Universitdt Heidelberg, Physikalisches Institut, Philosophenweg 12, D-69120 Heidelberg, Germany 

2 Now at Temple University, Philadelphia, Pennsylvania 19122, USA 
3 Also with Universita di Perugia, Dipartimento di Fisica, Perugia, Italy 



3 



F. U. Bernlochner, H. M. Lacker, T. Lueck, A. Volk 
Humboldt- Universitdt zu Berlin, Institut fur Physik, Newtonstr. 15, D- 12489 Berlin, Germany 

P. D. Dauncey, M. Tibbetts 
Imperial College London, London, SW7 2AZ, United Kingdom 

P. K. Behera, M. J. Charles, U. Mallik 
University of Iowa, Iowa City, Iowa 52242, USA 

J. Cochran, H. B. Crawley, L. Dong, V. Eyges, W. T. Meyer, S. Prell, E. I. Rosenberg, A. E. Rubin 
Iowa State University, Ames, Iowa 50011-3160, USA 

Y. Y. Gao, A. V. Gritsan, Z. J. Guo 
Johns Hopkins University, Baltimore, Maryland 21218, USA 

N. Arnaud, A. D'Orazio, M. Davier, D. Derkach, J. Firmino da Costa, G. Grosdidier, F. Le Diberder, 
V. Lepeltier, A. M. Lutz, B. Malaescu, P. Roudeau, M. H. Schune, J. Serrano, V. Sordini^ A. Stocchi, 

G. Wormscr 

Laboratoire de V Accelerateur Lineaire, IN2P3/CNRS et Universite Paris-Sud 11, Centre Scientifique 

d'Orsay, B. P. 34, F-91898 Orsay Cedex, France 

D. J. Lange, D. M. Wright 
Lawrence Livermore National Laboratory, Livermore, California 94550, USA 

I. Bingham, J. P. Burke, C. A. Chavez, J. R. Fry, E. Gabathuler, R. Garnet, D. E. Hutchcroft, D. J. Payne, 

C. Touramanis 

University of Liverpool, Liverpool L69 7ZE, United Kingdom 

A. J. Bevan, C. K. Clarke, F. Di Lodovico, R. Sacco, M. Sigamani 
Queen Mary, University of London, London, El 4NS, United Kingdom 

G. Cowan, S. Paramesvaran, A. C. Wren 

University of London, Royal Holloway and Bedford New College, Egham, Surrey TW20 OEX, United 

Kingdom 

D. N. Brown, C. L. Davis 
University of Louisville, Louisville, Kentucky 40292, USA 

A. G. Denig, M. Fritsch, W. Gradl, A. Hafner 
Johannes Gutenberg- Universitdt Mainz, Institut fur Kernphysik, D- 55099 Mainz, Germany 

K. E. Alwyn, D. Bailey, R. J. Barlow, G. Jackson, G. D. Lafferty, T. J. West, J. I. Yi 
University of Manchester, Manchester M13 9PL, United Kingdom 

J. Anderson, C. Chen, A. Jawahery, D. A. Roberts, G. Simi, J. M. Tuggle 
University of Maryland, College Park, Maryland 20742, USA 

C. Dallapiccola, E. Salvati 
University of Massachusetts, Amherst, Massachusetts 01003, USA 
4 Also with Universita di Roma La Sapienza, 1-00185 Roma, Italy 



4 



R. Cowan, D. Dujmic, P. H. Fisher, S. W. Henderson, G. Sciolla, M. Spitznagel, R. K. Yamamoto, M. Zhao 

Massachusetts Institute of Technology, Laboratory for Nuclear Science, Cambridge, Massachusetts 02139, 

USA 

P. M. Patel, S. H. Robertson, M. Schram 
McGill University, Montreal, Quebec, Canada H3A 2T8 

P. Biassoni ab , A. Lazzaro ab , V. Lombardo Q , F. Palombo Qb , S. Stracka ab 
INFN Sezione di Milano a ; Dipartimento di Fisica, Universita di Milano b , 1-20133 Milano, Italy 

L. Cremaldi, R. Godang0 R. Kroeger, P. Sonnek, D. J. Summers, H. W. Zhao 
University of Mississippi, University, Mississippi 38677, USA 

X. Nguyen, M. Simard, P. Taras 
Universite de Montreal, Physique des Particules, Montreal, Quebec, Canada H3C 337 

H. Nicholson 

Mount Holyoke College, South Hadley, Massachusetts 01075, USA 

G. De Nardo Qb , L. Lista Q , D. Monorchio ab , G. Onorato ab , C. Sciacca ab 

INFN Sezione di Napoli a ; Dipartimento di Scienze Fisiche, Universita di Napoli Federico IF, 1-80126 

Napoli, Italy 

G. Raven, H. L. Snoek 

NIKHEF, National Institute for Nuclear Physics and High Energy Physics, NL-1009 DB Amsterdam, The 

Netherlands 

C. P. Jessop, K. J. Knoepfel, J. M. LoSecco, W. F. Wang 
University of Notre Dame, Notre Dame, Indiana 46556, USA 

L. A. Corwin, K. Honscheid, H. Kagan, R. Kass, J. P. Morris, A. M. Rahimi, S. J. Sekula 
Ohio State University, Columbus, Ohio 43210, USA 

N. L. Blount, J. Brau, R. Frey, O. Igonkina, J. A. Kolb, M. Lu, R. Rahmat, N. B. Sinev, D. Strom, 

J. Strube, E. Torrence 

University of Oregon, Eugene, Oregon 97403, USA 

G. Castelli Qb , N. Gagliardi Qb , M. Margoni ab , M. Morandin a , M. Posocco a , M. Rotondo a , F. Simonetto ab , 

R. Stroili Qb , C. Voci ab 

INFN Sezione di Padova a ; Dipartimento di Fisica, Universita di Padova b , 1-35131 Padova, Italy 

P. del Amo Sanchez, E. Ben-Haim, G. R. Bonneaud, H. Briand, J. Chauveau, O. Hamon, Ph. Leruste, 

G. Marchiori, J. Ocariz, A. Perez, J. Prendki, S. Sitt 

Laboratoire de Physique Nucleaire et de Hautes Energies, IN2P3/CNRS, Universite Pierre et Marie 
Curie-Paris6, Universite Denis Diderot-P aris7 , F-75252 Paris, France 

L. Gladney 

University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA 
5 Now at University of South Alabama, Mobile, Alabama 36688, USA 



5 



M. Biasini afc , E. Manoni ah 
INFN Sezione di Perugia 0, ; Dipartimento di Fisica, Universitd di Perugia , 1-06100 Perugia, Italy 

C. Angelini ab , G. Batignani afc , S. Bettarini ab , G. Calderini ab @ M. Carpinelli ab Q A. Cervelli ab , F. Forti ab , 
M. A. Giorgi ab , A. Lusiani QC , M. Morganti ab , N. Neri ab , E. Paoloni ab , G. Rizzo ab , J. J. Walsh a 

INFN Sezione di Pisa a ; Dipartimento di Fisica, Universitd di Pisa b ; Scuola Normale Superiore di Pisa c , 

1-56127 Pisa, Italy 

D. Lopes Pegna, C. Lu, J. Olsen, A. J. S. Smith, A. V. Telnov 
Princeton University, Princeton, New Jersey 08544, USA 

F. Anulli a , E. Baracchini ab , G. Cavoto a , R. Faccini ab , F. Ferrarotto a , F. Ferroni ab , M. Gaspero ab , 
P. D. Jackson a , L. Li Gioi Q , M. A. Mazzoni a , S. Morganti a , G. Piredda a , F. Renga ab , C. Voena a 

INFN Sezione di Roma° ; Dipartimento di Fisica, Universitd di Roma La Sapienza , 1-00185 Roma, Italy 

M. Ebcrt, T. Hartmann, H. Schroder, R. Waldi 
Universitdt Rostock, D- 18051 Rostock, Germany 

T. Adye, B. Franek, E. O. Olaiya, F. F. Wilson 
Rutherford Appleton Laboratory, Chilton, Didcot, Oxon, 0X11 OQX, United Kingdom 

S. Emery, L. Esteve, G. Hamel de Monchenault, W. Kozanecki, G. Vasseur, Ch. Yeche, M. Zito 
CEA, Irfu, SPP, Centre de Saclay, F-91191 Gif-sur-Yvette, France 

M. T. Allen, D. Aston, D. J. Bard, R. Bartoldus, J. F. Benitez, R. Cenci, J. P. Coleman, M. R. Convery, 
J. C. Dingfelder, J. Dorfan, G. P. Dubois-Felsmann, W. Dunwoodie, R. C. Field, M. Franco Sevilla, 
B. G. Fulsom, A. M. Gabareen, M. T. Graham, P. Grenier, C. Hast, W. R. Innes, J. Kaminski, 
M. H. Kelsey, H. Kim, P. Kim, M. L. Kocian, D. W. G. S. Leith, S. Li, B. Lindquist, S. Luitz, V. Luth, 
H. L. Lynch, D. B. MacFarlane, H. Marsiske, R. Messnerd D. R. Midler, H. Neal, S. Nelson, 
C. P. O'Grady, I. Ofte, M. Perl, B. N. Ratcliff, A. Roodman, A. A. Salnikov, R. H. Schindler, 
J. Schwicning, A. Snyder, D. Su, M. K. Sullivan, K. Suzuki, S. K. Swain, J. M. Thompson, J. Va'vra, 
A. P. Wagner, M. Weaver, C. A. West, W. J. Wisniewski, M. Wittgen, D. H. Wright, H. W. Wulsin, 

A. K. Yarritu, C. C. Young, V. Ziegler 

SLAC National Accelerator Laboratory, Stanford, California 94309 USA 

X. R. Chen, H. Liu, W. Park, M. V. Purohit, R. M. White, J. R. Wilson 
University of South Carolina, Columbia, South Carolina 29208, USA 

M. Bellis, P. R. Burchat, A. J. Edwards, T. S. Miyashita 
Stanford University, Stanford, California 94305-4060, USA 

S. Ahmed, M. S. Alam, J. A. Ernst, B. Pan, M. A. Saeed, S. B. Zain 
State University of New York, Albany, New York 12222, USA 

A. Soffer 

Tel Aviv University, School of Physics and Astronomy, Tel Aviv, 69918, Israel 

8 Also with Laboratoire de Physique Nucleaire et de Hautes Energies, IN2P3/CNRS, Universite Pierre et Marie 
Curie-Paris6, Universite Denis Diderot-Paris7, F- 75252 Paris, France 
7 Also with Universita di Sassari, Sassari, Italy 
8 Deceased 



6 



S. M. Spanier, B. J. Wogsland 
University of Tennessee, Knoxville, Tennessee 37996, USA 

R. Eckmann, J. L. Ritchie, A. M. Ruland, C. J. Schilling, R. F. Schwitters, B. C. Wray 
University of Texas at Austin, Austin, Texas 18112, USA 

B. W. Drummond, J. M. Izen, X. C. Lou 
University of Texas at Dallas, Richardson, Texas 15083, USA 

F. Bianchi afc , D. Gamba ab , M. Pelliccioni afc 

INFN Sezione di Torino a ; Dipartimento di Fisica Sperimentale, Universita di Torino b , 1-10125 Torino, 

Italy 

M. Bombcn ab , L. Bosisio ab , C. Cartaro ab , G. Delia Ricca af> , L. Lanccri ab , L. Vitalc ab 
INFN Sezione di Trieste"", Dipartimento di Fisica, Universita di Trieste b , 1-34121 Trieste, Italy 

V. Azzolini, N. Lopez-March, F. Martinez- Vidal, D. A. Milanes, A. Oyanguren 
IFIC, Universitat de Valencia- CSIC, E-46071 Valencia, Spain 

J. Albert, Sw. Banerjee, B. Bhuyan, H. H. F. Choi, K. Hamano, G. J. King, R. Kowalewski, 
M. J. Lewczuk, I. M. Nugent, J. M. Roncy, R. J. Sobie 

University of Victoria, Victoria, British Columbia, Canada V8W 3P6 

T. J. Gershon, P. F. Harrison, J. Ilic, T. E. Latham, G. B. Mohanty, E. M. T. Puccio 
Department of Physics, University of Warwick, Coventry CV4 1AL, United Kingdom 

H. R. Band, X. Chen, S. Dasu, K. T. Flood, Y. Pan, R. Prcpost, C. O. Vuosalo, S. L. Wu 
University of Wisconsin, Madison, Wisconsin 53106, USA 



7 



1 Introduction 



Recent cosmic ray measurements of the electron and positron flux from ATIC[T], FERMI [2 j, and 
PAMELA [3] have spectra which are not well described by galactic cosmic ray models such as 
GALPROP[4]. For instance, PAMELA shows an increase in the positron/electron fraction with 
increasing energy. No corresponding increase in the antiproton spectrum is observed. There have 
been two main approaches attempting to explain these features: astrophysical sources (particularly 
from undetected, nearby pulsars) [5j and annihilating or decaying dark matter. 

Arkani-Hamed et al.\£>\ have introduced a class of theories containing a new "dark force" and 
a light, hidden sector. In this model, the ATIC and PAMELA signals are due to dark matter 
particles with mass ~ 400 — 800 GeV/c 2 annihilating into the gauge boson force carrier with mass 
~ lGeV/c 2 , which they dub the eft, which subsequently decays to Standard Model particles. If 
the 4> mass is below twice the proton mass, decays to pp are kinematically forbidden allowing only 
decays to states like e + e~, , and tttt. If the dark force is non-Abelian, this theory can also 

accommodate the 511 keV signal found by the INTEGRAL satellite [7] and the DAMA modulation 
data [8]. 

The dark sector couples to the Standard Model through kinetic mixing with the photon. Thus 
low-energy /high luminosity e + e~ experiments like BABAR are in excellent position to probe these 
theories. Recent papers by Batell et al. [9j and Essig et al. [TO] have discussed the prospects 
for finding evidence for the dark sector at the B-Factories in the Abelian and non-Abelian cases, 
respectively. In the Abelian case, the signatures would be e + e~ — > j(p ~* or e + e~ — > 

<j)h' — ► 3(/ + /~) (where h' is a "dark Higgs"). There are actually two non-Abelian scenarios: the 
Higgsed case and the confined case ("dark QCD"). In the Higgsed case there are at least three 
dark particles in play: A' which mixes with the photon, another gauge boson W', and the dark 
Higgs b! . In this regime, signatures are e + e~ — > WW — * l + l~l + l~ (via a virtual A') and e + e~ — > 
yA'(—> WW) — > , plus "Higgs'-strahlung" processes which may lead to missing energy. 

Finally, the confined case could lead to a proliferation of "dark mesons" , whose lowest mass states 
decay to leptons. Depending on the scenario and the coupling between the Standard Model and 
dark sectors, cross sections could be as large as a few femtobarns at BABAR which would translate 
to hundreds of events observed in the detector. 

In this note we describe a search for the W' in the reaction e + e~ — > WW — ► in 
exclusive 4-lepton final states, where we require that the four leptons carry the full center of mass 
energy and that the two dilepton pairs have the same invariant mass. 

2 The BABAR Detector and Dataset 

The data used in this analysis were collected with the BABAR detector at the PEP-II asymmetric 
energy e + e~ storage rings between 1999 and 2008 and correspond to an integrated luminosity of 
536 fb . This data was mostly at the T(4S) peak but it also includes collisions at the T(2S) and 
T(3S) as well as off-resonant data. 

To study signal efficiency and resolution, e + e~ — > WW — ► l + l~l + l~ Monte Carlo (MC) 
samples were generated (where l=e or fj,) for different values of W mass using the MadGraph 
event generator [TT|. There were 10 4 events generated at each mass value of: 0.3, 0.4, 0.5, 0.7, 1.0, 
1.5, 2.0, 3.0, 4.0, and 5.0 GeV/c 2 . To study backgrounds, we have inspected BB (~3x luminosity), 
uds, cc, and rr MC samples (each ~lx luminosity). In addition we created 4-lepton QED samples 
using the diag36 event generator [12J. 
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A detailed description of the BABAR detector is given in [13] . Charged-particle trajectories are 
measured by a five-layer, double-sided silicon vertex tracker (SVT) and a 40-layer drift chamber 
(DCH) coaxial with a 1.5 T magnetic field. Charged-particle identification is achieved by combining 
the information from a ring-imaging Cherenkov device (DIRC) with the ionization energy loss 
(dE/dx) measurements from the DCH and SVT. Photons are detected in a CsI(Tl) electromagnetic 
calorimeter (EMC) inside the coil. Muon candidates are identified in the instrumented flux return 
(IFR) of the superconducting solenoid. We use GEANT4-based [H] software to simulate the 
detector response and account for the varying beam and environmental conditions. 

3 Event Selection 

We search for the exclusive pair production of a narrow resonance, consistent with the detector 
resolution, decaying to leptons and with a mass in the range between 240MeV/c 2 to y/s/2. The 
signature is 4 leptons with zero total charge carrying the full beam momentum where the two 
dilepton invariant masses are equal. This topology, particularly the equal invariant masses, is quite 
unique and the only backgrounds are from 4-lepton QED processes. The full selection criteria are 
described below. We used 10% of the data as a test (blind) sample to choose our selection and 
signal extraction procedures before looking at the full dataset. 
We begin by selecting events with: 

• 4 charged tracks 

• two leptons with pcm > 1.5GeV/c 

• sum of the absolute value of momentum of all tracks > 6GeV/c 2 or the total visible energy 
(lab)> 8GeV/c 2 

We reconstruct 4-lepton candidates from combinations of two W — * l + l~ candidates. The 
lepton candidates are chosen by their signitures in the EMC and IFR. The W candidates are 
formed from e + e~ or pairs. We then select events which satisfy the following criteria: 

. [iV e ,Ay = [4,0], [2, 2], or [0,4] 

• M Alepton > 10GeV/c 2 

• the helicity angle of a lepton pair, defined as the angle between the positive lepton and the 
lepton-pair flight direction, is required to be \cos(9h)\ < 0.95 for each pair 

• to reduce background from photon converstions, we require the flight significance, defined as 
the W candidates decay length from the interaction point divided by the error, is < Act for 
each pair 

• to reduce background from radiative Bhabha events, we require the angle between the decay 
planes of the lepton pairs, <f>DPN > 0.2 

The 4-lepton candidate is then fit constraining the four-momentum to the total beam momentum 
and the vertex to the interaction point. 

At this point, we can exploit the fact that both dilepton pairs for our signal events have the 
same invariant mass. The 2-dimensional distributions of dilepton masses for each final state after 
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Figure 1: The dilepton invariant mass distributions from data for (left to right) e + e e + e , 
e + e~u, + u,~, and after all other cuts. The solid lines denotes m-i = mi- 

all of the above cuts for the blinding sample is shown in Figure HJ We define the transformed 
masses: 

m = (mi + m 2) /2 (1) 

Am = | mi — m,2\ (2) 

where m\ and ni2 are the dilepton invariant masses. The distribution of events for these variables 
is shown in Figure [2j We impose a cut on Am (shown as the solid line in Figure [2|) of Am < 
0.25GeV/c 2 for m < 1.0GeV/c 2 and Am < 0.50GeV/c 2 for m > 1.0GeV/c 2 . Because of the 
threshold effects in n, + n,~ n, + yT , we tighten the Am cut in a linear fashion below m < 4 x M(n). 

In the case of the e + e~e + e~ and u, + ur u, + n,~ final state, there are two possible l + l~ pair 
combinations. If both pairings pass all cuts, the pair with the smallest value of Am is used. For 
data, we see two pairings passing all cuts except the Am cut for 25% of e + e~e + e~ events and for 
44% of the events. Table CD shows the progressive and total efficiencies for the three 

different final states of WW — > (assuming the mass of the W is lGeV/c 2 ) as well as 

the progressive efficiency for the data. As shown in the table, the loose cut on Am is extremely 
powerful at reducing the background while not affecting the signal efficiency. After all selection, 
there are 28303 events remaining in our data sample; of these 16531 are e + e~e + e~ events, 9592 
are e + e~/x + /i~ events, and 2180 are events. 



4 Signal Extraction 

Our aim is to perform a search for a narrow peak in the m range from 240MeV/c 2 up to \fsj2. 
After the selection described in the previous section, the expected backgrounds are quite low and 
we have decided to perform a cut-and-count analysis in bins of to, using the the Am variable to 
define the signal and background regions. The number of observed signal events in a to bin is then: 

N s ig — -^signal region -^bkg region X — ■ (3) 

-^background 

where A signa i (A^ ac k groun( i) is the area of the signal (background) Am region. 

In this section, we will discuss the signal efficiency, Am shapes (including the definition of signal 
and background regions) and background rates as a function of m and the method we plan to use 
in extracting the signal yields and setting limits. 
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Figure 2: The transformed mass distributions, Am vs m, from data for (left to right) e + e 
e + e~ and n + fi~fi + n~ after all other cuts. The solid lines denotes the Am cut value. 
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4.1 Efficiency and Am resolution dependence on the W mass 

The efficiency for different generated values of the W' mass is shown in Figure [3l The efficiency 
decreases from ~ 45% at 1 GeV to 25 — 30% at high masses depending on the decay mode. There is 
a dip in efficiency for mass pairs around 500 MeV/ c 2 which is due to the opening angle of the lepton 
pair at this mass coinciding with the bending angle at the EMC, precluding us from identifying the 
two particles for a fraction of the events. The Am resolution also varies significantly as a function 
of W mass. Figure H] shows the distributions of Am for four different mass values. The resolution 
of Am increases with increasing W' mass. Since the background Am distribution is basically flat 
and roughly constant in m (see Section T4.2p . the effect is to reduce the sensitivity at higher masses. 

Figure [5] shows the values of the Am cut which retains 90% of the signal as a function of m. 
We use this cut value to define the signal (Am < cutVal) and background (dm > cutVal) regions 
for the cut-and-count signal extraction. Recall that the maximum value of Am is 0.25(0.5) GeV/c 2 
for m < (>)1.0 GeV/c 2 . The solid line is the result of a 4th-order polynomial fit which we use to 
extrapolate between m points. 

4.2 Background composition 

While we ultimately use the Am sidebands to determine our background level, we have also used 
MC to study the composition of the background. In generic qq, B°73° , B + B~ , and t + t~ samples 
we find only a single event passing the cuts (a qq event in the 4-electron final state). From this we 
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Table 1: Selection efficiencies relative to the previous cut with binomial errors for the three signal 
decay modes assuming M(W') = 1 GeV/c 2 and for onpeak data. 



Cuts Relative Efficiencies (%) 

£ W'W'^>Ae £ W'W'^2e2/i e WW^4/i £ data 



Preselection 


61.4 ± 


0.5 


68.2 ± 


0.5 


73.4 ± 


0.4 








N(tracks)=4 


93.3 ± 


0.3 


95.5 ± 


0.3 


97.1 ± 


0.2 


83.7 


± 


0.0 


N(leptons) 


99.9 ± 


0.1 


100.0 ± 


0.0 


100.0 ± 


0.0 


98.7 


± 


0.0 


M (41) > 10 GeV 


87.2 ± 


0.5 


93.2 ± 


0.3 


98.1 ± 


0.2 


66.2 


± 


0.0 


\cos(6 H )\ < 0.95 


99.9 ± 


0.1 


99.6 ± 


0.1 


98.8 ± 


0.1 


18.2 


± 


0.0 


fit. sig.< 4<T 


97.2 ± 


0.2 


96.5 ± 


0.4 


99.1 ± 


0.1 


54.6 


± 


0.1 


fj, PID 


100.0 i 


0.0 


81.9 ± 


0.6 


70.3 ± 


0.6 


74.2 


± 


0.1 


<t>DPN > 0.2 


92.4 ± 


0.4 


93.2 ± 


0.4 


93.6 ± 


0.4 


37.0 


± 


0.1 


A(m) 


98.3 ± 


0.2 


99.8 ± 


0.1 


99.6 ± 


0.1 


3.8 


± 


0.1 


Total Efficiency 


43.7 ± 


0.5 


44.5 ± 


0.5 


44.8 ± 


0.5 
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Figure 3: The signal efficiency versus W' mass for (left to right) W'W — > e + e e + e 
e + e~ / u + /i~, and WW — > fi + fi + after all cuts. 



W'W 



conclude that our background is dominated by QED processes. 

We have generated e + e~/x + /i~ and /j + fi~^ + fi~ sampled using the diag36 generator and com- 
pared the MC to our selected dataset. We find good agreement both in the scale and shape between 
data and the four-lepton QED MC. From the MC, we expect to observe 16241 ± 250 fi + fi~fi + fi~ 
events in the full dataset while we observe 15666 ± 125 (statistical errors only). For e + e~/x + ^ _ we 
expect 219927 ± 3450 and observe 185499 ± 431 events. 

The background distributions in Am and m, after all selection, are shown in Figure [6j The 
background Am distributions were fit with a line in different slices of m, the slopes of which are 
plotted for the three modes in Figure [71 and the slopes are consistent with 0. When extracting the 
signal yields, we assume a uniform background distribution, and take into account the uncertainties 
in the slope as a systematic error. We use the full dataset for the above plots; any signal present 
would be completely washed out when projected onto the Am or m axis. 

9 Due to the enormous e + e~e + e~ QED cross-section, this mode is difficult to generate efficiently. 



12 



BaBar 
preliminary 



■ 500 MeV/c 2 
1 1000 MeV/c 2 
3000 MeV/c 2 
1 5000 MeV/c 2 



0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 

A m (GeWc 2 ) 



BaBar 
preliminary 



■ 500 MeV/c 2 
1 1000 MeV/c 2 
3000 MeV/c 2 
1 5000 MeV/c 2 



0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 

A m (GeWc 2 ) 



BaBar 
preliminary 



■ 500 MeV/c 2 
1 1000 MeV/c 2 
3000 MeV/c 2 
1 5000 MeV/c 2 



°0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 "6*08 o'.09"'0.1 

A m (GeWc 2 ) 



Figure 4: The Am distributions for four different W mass values (left to right) WW 
WW —* e + e~/i + /i~, and WW —* after all cuts. 
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Figure 5: The values of the cut on Am keeping 90% of signal events as a function of W mass for 
(left to right) WW -> e+e'e+e', WW e + e"^ + /j-, and WW -> The line is a 

fit to a fourth order polynomial. This cut defines our signal and background region. 



4.3 Signal extraction and limit setting 

In this analysis, our aim is to obtain a limit (or observe a signal) for e + e~ — > W'W' as function 
of the presumed W mass. To this end, search for a signal in steps of the average dilepton mass 
m. We have chosen the m bin size to be 20MeV/c 2 , which is a large enough range to fully contain 
any signal. Figure [8] shows the RMS of m at the different mass points. We scan m in steps of 
10MeV/c 2 , half the bin size, so that at least one bin will fully contain the signal. Thus, in the m 
range from 0.24 — 5.3 GeV/c 2 , there are 507 total bins. We define the signal and background regions 
in Am by cutting at a value of Am so that the signal is 90% efficient, as discussed above. 

With this framework, the number of background events in a given m bin is quite small. Except 
at low m, the expected number of background events in the entire Am range is typically below 100 
events in a m bin, particularly for the [i + fi~fi + fi~ mode where it is below 5 events. Thus there will 
be relatively large fluctuations in the background due to Poisson statistics and the limit setting 
procedure must take this into account. We use a profile likelihood technique [15] to set limits in the 
presence of nuisance parameters, such as the expected background yield. Using this technique, we 
obtain a confidence level (CL) for the presence of signal defined as: 

CL = Prob (-2 log(£ s=0 ) - 21og(£ ma ,)) (4) 

where £ s =o is the value of the likelihood at signal events and C ma x is the maximum value of the 
likelihood. 

Since in our dataset we will have 507 correlated measurements (204 independent measurements), 
each at a different m, we need to determine a criteria for a signal observation. Simply asking whether 
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Figure 6: The background (top) Am and (bottom) m distributions for (left-to-right) e + e - e + e~, 
e + e~ fi + fi~ , and fi + [i~ from the full dataset. For the Am plots, we have required m > lGeV 
to that all events have the same upper Am value. The effect of the Am cut increasing from 
0.25GeV/c 2 to 0.5GeV/c 2 at m = 1.0GeV/c 2 can be seen in the m plots. 



an individual bin has an observed yield in it > 3a above is not enough since the probability to 
observe at least 1 > 3<r fluctuation in one of the m bins is 0.3 (as determined from the simulation 
described below). We need to redefine the Xa levels for the new question "What is the chance 
that I see a background fluctuation above Xa in our 507 correlated trials?" . We have done this by 
generating many simulated datasets (toys) with the expected m and Am background distributions 
with signal and plotting the highest value of the signal confidence level observed over that dataset, 
which we call CL max . The results of these simulations are shown in Figure [9l plotting the more 
convenient variable —ln(l — CL max ). As a reference, the distribution of values —ln(l — CL) from 
a single bin (i.e. not the largest value in an m scan) is shown in Figure PTOl Table [2] shows the 
values of — ln(l — CL max ) that correspond to l-4a fluctuations of the background (also displayed 
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Figure 7: The slope of the background Am distributions for (left-to-right) e + e e + e 



and as a function of m. The mean values of the slopes are: 

and -0.02 ±0.19. 



-0.11 ±0.06, 0.07 ±0.08, 
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Figure 8: The m RMS versus W mass for (left to right) W'W 
and WW -» fi + fi- fi + fi- . 
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on the plot). Although the background levels are different, the values are consistent between the 
three modes. 

Additionally, we calculate the combined max confidence level, defined as: 



(1 — CL max ^c) — (1 — CL maXt 4e)(l — CL 



max,2e2fi 



)(l-CL 



(5) 



whose distribution for background-only toys is shown in the bottom right plot of Figure [9l If 
lepton universality holds, this limit is potentially more sensitive than the individual confidence 
levels and allows us to catch a signal that is not significant in any single final state. Our criteria to 
claim evidence of a signal is to observe the largest value of — ln(l — CL max ) in any of e + e~e + e~, 
e + e _ / u + /i~, fi + fj,~ or in the combined confidence level that is greater than the 3a values given 
in Table H 



Table 2: Values of the 1-,2-,3-,4-<t limits for — ln(l — CL max ) in the three final states. 



Signif. P(CL max ) < X -Zn(l - CL max ) 







e + e e + e 


e + e 




Combined 


la 


0.84135 


7.2 


7.1 


7.0 


10.2 


2a 


0.97725 


9.3 


9.2 


9.0 


12.4 


3(7 


0.99865 


12.2 


12.1 


11.6 


15.8 


la 


0.99997 


16.3 


16.1 


14.5 


19.2 



5 Systematic Errors 

There are two types of systematic errors in this analysis: systematics that effect both the yield and 
cross-section upper limits (e.g. errors due to uncertainties in the background shape) and systematics 
that just affect the cross-section (e.g. tracking efficiency errors). The second type of error does not 
effect the signal significance. Table [3] summarizes the values the systematic errors for the different 
sources described below. 

• Am background shape: We assume that the background is uniform in Am and with our 
limited MC statistics but we have no a priori reason to expect this. While the background 
Am does look quite flat and does not appear to depend on m, see Figure EJ, we still need 
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Figure 9 

C L max ) 



; The distribution of —ln(\ — CL max ) from toy with the arrow showing the value of —ln{\ — 
observed in data. The plots are (left to right, top to bottom) e + e~e + e~, e + e _ /i + /i _ , 
and the three modes combined. 



to account for uncertainties. Consequently, we estimate the Am background shape from the 
data itself. 

In order to estimate the size of this uncertainty, we have generated toy m scans (background 
only) with a slope and calculated the signal yield assuming a slope of 0. The mean Am slopes 
are given in the caption of Figure For this study, we shift the mean value of the slope (B m ) 
by: 

— for B m < 0, we assign the slope to be B m — a 

— for B m > 0, we assign the slope to be —a 

where a m is the error on the mean. We only use the negative slope values because we are 
primarily interested in how this biases us toward more signal. The results of this study are 
shown in Figure Qj] as the observed signal yield bias vs m for the three modes. The bias 
depends on m because both the number of background events in the full Am region and 
because the Am signal/background region definitions depend on m. 

We incorporate this bias into a systematic error on the cross section by converting the bias 
in the number of events into a cross section in m bins. The error is largest for the e + e~e + e~ 
mode where at high m is as large as ~ 5 ab; for the other two modes this error is generally 
< 1 ab. 
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Figure 10: The distribution of values of — ln(l — CL) from (error bars) data and (solid histogram) 
toy for (left to right) e + e~e + e~, e + e~fi + fi~, and fi + \i~ /i + \T~ . 
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Figure 11: The positive signal yield bias due to the uncertainty in the background Am slope as a 
function of m for (left-to-right) e + e~e + e~, e + e~ fj, + fj,~ , and 



• Am signal shape: We use MC at select mass values to interpolate the 90% efficiency Am 
cut value to cover all masses. The interpolation is done with a polynomial and we vary the 
parameters of the polynomial within their errors to get the error in the Am cut value. This 
is then translated into an efficiency error. The magnitude of this error is ~ 1% and depends 
slightly on m. 

• interpolation of total efficiency: We use MC at select mass values to interpolate the 
total efficiency to all masses. The interpolation is done by interpolating the efficiency linearly 
between the MC mass points. We propagate the errors in the efficiency points due to MC 
statistics through the interpolation. In addition, we take the difference between a linear and 
quadratic interpolation and assign the difference, added in quadrature with the statistical 
error, as the systematic. The magnitude of this error is ~ 3% and depends slightly on m. 

• particle ID: we assign a 1% error per electron and 2% error per muon on the cross sections 
to account for the systematic error in the PID efficiency. This is the dominant systematic 
error. 

• tracking efficiency: we assign 0.21% error per track on the cross sections to account for 
the systematic error in the charged track reconstruction efficiency. 

• luminosity: we assign a 1.1% error on the cross sections due to the uncertaintly in the total 
luminosity. 

We add these sources of systematic error in quadrature and scale the statistical 90% upper limit 
by the fractional systematic error to obtain the final upper limit. This error depends slightly on 
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m but is around 5% for e + e e + e (up to 10% for high m), 6.5% for e + e , and 8.2% for 

n+p-n+ijr. 

Table 3: Sources of systematic uncertainties and their contributions. 
Source e + e~e + e~ e + e~fj, + /j,~ \x~ /x + \i~ 

Am bkg shape 0.4-5.5 ab 0.1-0.7 ab 0.1-0.3 ab 



Am signal efficiency 1% 1% 1% 

total signal efficiency 3% 3% 3% 

particle ID 4% 6% 8% 

tracking efficiency 0.8% 0.8% 0.8% 

luminosity 1.1% 1.1% 1.1% 



6 Results and Conclusions 

The spectra for the entire dataset (including the 10% test sample) show no significant signal in 
any of e + e~e + e~, e + e~ / u + //~, /i + / u~ / u + /^ _ final states, or the combination of the three. The 
summary of results is shown in Tabled! The distribution of observed signal events, after background 
subtraction, for all bins in m is shown in Figure [12] The values of — ln(l — CL) versus m are shown 
in Figure [13] and show no bins above the 3a value, shown on the plots. The raw distribution of 
— ln(l — CL) compared to toy simulations with only background is shown in Figure fTUl and is in 
good agreement. The plots in Figure [9] compare the values of tlie — Iti{\. — GL rnQjX ) observed in data 
with the distribution found in toy simulation. 

Table 4: Summary of the — ln(l — CLmax) observed in data. 





— ln{\ — CL max ) 


m max (GeV) 


e + e e + e 


5.88 


5.27 


e + e~ \i r ii~ 


6.26 


1.44 




6.94 


2.23 


Combined 


7.15 


1.66 



Correcting for efficiency (Figure El using linear interpolation between points and including the 
90% cut on Am) and scaling by the luminosity, we obtain a 90% upper limit for the cross section as 
shown in Figures [14] and [15] The points in these plots are the upper limit for each bin in m while 
the solid lines are the averages of the upper limits in the m region shown. We set upper limits of 
a(e + e~ -> WW -► e+e'e+e") < (15 - 70) ab, a(e + e~ -» WW -» e + e~^ + ^-) < (15 - 40) ab, 
and cr(e + e~ — > WW — ► fi + [T~ fi + ^) < (11 — 17) ab depending on W' mass (taking the ranges 
from the averaged limits). 

Assuming lepton universality (BR(W — > e + e~) = BR(W — > /i + /u~)), we combine the three 
modes to obtain upper limits for the reaction e + e~ — > WW — * . We obtain this limit 

by combining the individual profile likelihood functions for the three decay modes as a function of 
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Figure 12: The number of signal events after background subtraction versus m for (left to right) 
e + e~e + e~, e + e~fi + fi~, and fi + /i~ fi + /i~ . The band structure evident in the fi + fi~ fi + fi~ plot is due 
to the very low number of events in this mode. 
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Figure 13: The value of — ln(l — CL) versus m for (left to right) e + e e + e , e + e , and 



e + e — > WW — > Z + Z cross section. The combined upper limit is shown in Figure [T5l we set 

upper limits for a(e + e~ -> W'W -> Z+n'+Z' - < (25 - 60) ab. 

From the combined upper limit, we derive limits on the possible couplings between the Standard 
Model and dark sectors. The cross section for e + e~ — ► W'W has been calculated by Essig et aLflOj. 
For a dark photon .A' mass less than the center of mass energy, E cm , the cross section is given by: 

, 47r e 2 ana / 4m?,,, / 2m?m \ , , 

CT(e + e - ^ w'w')^ = N^-f^h\ 1 - -trF 1 + (6) 

" ^cm y cm \ *- J cm J 

while for an A' mass larger than £" cm the cross section is: 

,(«+.- - ^ W = Nc f^^l-^ + (7) 

where iV c is the number of colors in the dark sector, e is the mixing parameter between the SM and 
the dark sector, and an is the dark sector coupling constant. Figure [TBI shows the upper limits we 
obtain on e 2 ar> assuming low A' mass, or on e ° D assuming large A' mass. For most of the mass 

m A , 

range, we exclude values of e 2 aij above 2 x 10~ 10 in the low A' mass scenario or values of £ 4 

m A> 

above 2 x 10 -14 in the high A' mass scenario. In the model of Ref |10j . these limits exclude the 
preferred parameter region for A' masses above 1.0GeV/c 2 . 
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Figure 14: The cross section 90% upper limit versus m for (top to bottom) e + e~ — > WW — > 
e + e~e + e~ and e + e~ — > W'W' — ► e + e~/i + /i~. The points are the upper limit for each m bin while 
the lines are the average of the limits over many bins. 
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Figure 15: The cross section 90% upper limit versus rn for (top to bottom) e + e~ — > WW' — > 
fi + fi~fi + /i~ and the combined e + e~ — > WW — > assuming lepton universality. The 

points are the upper limit for each m bin while the lines are the average of the limits over many 
bins. The band structure evident in the plot is due to the very low number of events 

in this mode. 
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Figure 16: The 90% upper limit on e 2 ao (left axis) or Mr 2 - (right axis) versus m(W). The points 



A' 



are the upper limit for each m(W') bin while the lines are the average of the limits over many bins. 
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